The temperature dependent behaviour of surface states in ferromagnetic semiconductors 
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We present a model calculation for the temperature dependent 
behaviour of surface states on a ferromagnetic local-moment film. 
The film is described within the s-f model featuring local magnetic 
moments being exchange coupled to the itinerant conduction elec- 
trons. The surface states are generated by modifying the hopping 
in the vicinity of the surface of the film. In the calculation for the 
temperature dependent behaviour of the surface states we are able to 
reproduce both Stoner-like and spin-mixing behaviour in agreement 
with recent (inverse) photoemission data on the temperature depen- 
dent behaviour of a Gd(OOOl) surface state. 



In the recent past many theoretical and experimental re- 
search works have been focussed on the intriguing properties 
of rare earth metals and their compounds. On the exeperi- 
mental side, this interest was aroused after Weller et at. re- 
ported on the existence of the magnetically ordered Gd(OOOl) 
surface at temperatures where the bulk Gadolinium is para- 
magnetic [|l]]. Since then, a variety of different experimen- 
tal techniques has been applied to the problem by different 
groups yielding values of the surface Curie temperature en- 
hancement, ATc = Tc(surface) - Tc(bulk), between 17K 
and 60K Contrary to the groups cited above Donath 

et at. using spin-resolved photoemission did not find any in- 
dication for an enhanced Curie temperature at the Gd(OOOl) 
surface [|]|, fueling the controversial discussion. 

Concerning the interplay between electronic structure and 
exceptional magnetic properties at the Gadolinium surface a 
Gd(OOOl) surface state is believed to play a crucial role 
and its temperature dependent behaviour has been discussed 
intensely [^^[jis]]. Recently, the investigation of the coiTe- 
lation between strain induced alteration of the surface elec- 
tronic structure and enhanced magnetization in Gd films has 
been addressed by a number of experimental works [13- 1^] . 
A thorough account on the surface magnetism and the surface 
electronic structure of the lanthanides has been given by Dow- 
ben ef a/. [|l7]|. 

Rare-earth materials are so-called local-moment systems, 
i.e. the magnetic moment stems from the partially filled 4f- 
shell of the rare-earth atom being strictly localized at the ion 
site. Thus the magnetic properties of these materials are de- 
termined by the localized magnetic moments. On the other 
hand, the electronic properties like electrical conductivity are 
borne by itinerant electrons in rather broad conduction bands, 
e.g. 6s, 5d for Gd. Many of the characteristics of local- 



moment systems can be attributed to a correlation between 
the localized moments and the itinerant conduction electrons. 
For this situation the s-f model has been proven to be an ad- 
equate description. In this model, the correlation between lo- 
calized moments and conduction electrons is represented by 
an intraatomic exchange interaction. 

In what follows we consider a film consisting of n equiva- 
lent parallel layers. The lattice sites within the film are indi- 
cated by a greek letter a, /3, 7, . . ., denoting the layer, and by a 
latin letter i, j, k, . . ., numbering the sites within a given layer. 
Each layer posseses two-dimensional translational symmetry, 
so for any site dependent operator Aia we have: 

The Hamiltonian for the i-/ model consists of three parts: 



7i = "Hs + "Hf + Ti-sf- 



The first. 
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describes the itinerant conduction electrons as s-electrons 
with (fiiaa) being the creation (annihilation) operator of 
an electron with the spin a at the lattice site RiQ. T'^^^ are the 
hopping integrals. 

The second part of the Hamiltonian represents the system 
of the localized /-moments and consists itself of two parts. 
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where the first is the well-known Heisenberg interaction. Here 
the Sia are the spin operators of the localized magnetic mo- 
ments, which are coupled by the exchange integrals, J'^f ■ The 
second contribution is a single-ion anisotropy term which ar- 
rises from the necessity of having a collective magnetic or- 
der at finite temperatures, T > [p^. This anisotropy has 
been assumed to be uniform within the film. The accord- 
ing anisotropy constant Dq is typically smaller by some or- 
ders of magnitude than the Heisenberg exchange integrals, 
« Jtf. 

In addition to the contribution of the s-electron system and 
the contribution of the localized /-moments we have a third 
term which accounts for an intraatomic interaction between 
the conduction electrons and the localized /-spins: 
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where J is the s-f exchange interaction and aia is the Pauh 
spin operator of the conduction electrons. In the case where 
J < the Hamiltonian (|l]) is that of the so called Kondo lat- 
tice. However, here we are interested in the case of positive 
s-f coupling ( J > 0) which applies to the materials we are 
interested in. Using the abbreviations, 

^ia = ^ia + l-^crSf^; = ±1, 

the s-/Hamiltonian can be written in the form: 

J • 



■ {zaS^^niaa + S^^cf^_^Cioia) ■ (5) 



The problem posed by the Hamiltonian can be solved 
by considering the retarded single-electron Green function 



E 







, (6) 



which is related to the spectral density S^^{E) and the local 
density of states paa{E) via the relations: 



S'^^{E) = -^lu,GZ{E + iQ+\ 



(7) 

(8) 
(9) 



Due to the translational symmetry of the films, the Fourier 
transformation (0) has to be performed within the layers of 
the film. Accordingly, N is the number of sites per layer, k 
is an in-plane wavevector from the first two-dimensional Bril- 
louin zone, and represents the in-plane part of the position 
vector, = + r^^. 

The many-body problem that arises with the Hamiltonian 
([l|) is far from being trivial and a full solution is lacking even 
for the case of the bulk. In a previous paper we have pre- 
sented an approximate treatment of the special case of a sin- 
gle electron in an otherwise empty conduction band [|l9||. This 
solution, which holds for arbitrary temperatures, is based on 
the special case of an empty conduction band interacting with 
a ferromagnetically saturated local-moment system (T = 0), 
which can be solved exactly for both, the bulk case and the 
film geometry [pO[-p^]. This exactly soluble limiting case 
gives the approximate solution for finite temperatures a cer- 
tain trustworthiness. 

In the following formulas, we briefly recafl the results of the 
calculations presented in | ]r^ , p3[ |. Due to the empty conduc- 
tion band, we are considering throughout the whole paper, the 



Hamiltonian ([l|) can be split into an electronic part, Tis+'Hs f, 
and a magnetic part, TC f, which can be solved separately Jl^]. 
For the electronic part. 



(10) 



we employ the single-electron Green function (^. The equa- 
tion of motion for this Green function G^^(£') can be for- 
mally solved by introducing the self-energy M"^ (E), 



J2MZJE)G':f^AE), (11) 



which contains all the information about correlation between 
the conduction band and the localized moments. With the help 
of ( pi] ) the equation of motion for the single-electron Green 
function simply becomes, after two-dimensional Fourier 
transform. 



(12) 



Here, I represents the {n x n) identity matrix and the matri- 
ces G]ia-{E), Tk, and Mko-(£') have as elements the layer- 
dependent functions G^^{E), T^^ , and M^^{E), respec- 
tively. The necessary computation of the self-energy M^^{E) 
involves the evaluation of higher Green functions originating 
from the equation of motion of the original single-electron 
Green function. For the sake of brevity we here omit the de- 
tails of the calculations. If the self-energy is assumed to be a 
local entity. 



M^!{E) = -"-S^^m^E), 
it can be shown to have the structure: 



J, 



N^iEY 



(13) 



(14) 



where the numerator and the denominator, respectively, have 
the structure: 

Z: = z,h^{Sl) + ^/f (. . .) + . .) (15a) 
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(15b) 



where the four functions fi2^ {■ ■ ■) themselves depend |^ on 
the self-energy, m'^^{E), and on the layer-dependent coeff- 
cients: 



1) 



A(2) 



{Sa S^) 



(16) 



^For the explicit form of Eqs. {MSl) see ni5 
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Taking into account Eqs. ( 12 )-(|l6|), we have now a closed sys- 
tem of equations, provided that the /-spin correlation functions 
appearing in Eqs. ( [T^ ) are known. 

These can be evaluted by considering the magnetic subsys- 
tem, according to the Hamiltonian Hf (Eq. (^). For all what 
concerns us in this paper it is only interesting that, employing 
an RPA-type decoupling, a solution of the magnetic subsys- 
tem can be found [|l8]| and that it gives us all the necessary 
layer-dependent /-spin correlation functions as a function of 
temperature from T = to the Curie temperature, T = Tc 
[p^]. Mediated by Eqs. ( 14)-(|l6|), the /-spin correlation func- 
tions contain the whole temperature dependence of the elec- 
tronic subsystem. 

To briefly recall the main results presented in our previous 
paper [|l9|]. Fig. |l] shows the density of states of a s.c.-(lOO) 
double layer for different s-f interactions and different tem- 
peratures. In the case of ferromagnetic saturation, T = 0, we 
see that for the spin-t electron the density of states of the free 
case J = (dotted line) is just rigidly shifted when the in- 
teraction is switched on. This is due to the impossibility of 
the spin-t electron to exchange its spin with the perfectly al- 
ligned local-moment system. For the spin-| electron in the 
case of small s-f exchange coupling, J > 0, a slight defor- 
mation of the free density of states sets in. For intermediate 
and strong couplings (J ^ 0.2), the density of states splits 
into two parts, corresponding to two different spin exchange 
processes between the spin-| electron and the localized /-spin 
system. The higher energetic part represents a polarization 
of the immediate spin neighbourhood of the electron due to 
the repeated emission and reabsorption of magnons. The cor- 
responding polaron-like quasiparticle is called the magnetic 
polaron. The low-energetic part of the spectrum is a scatter- 
ing band which can be explained by the simple emission of a 
magnon by the spin-^ electron without reabsorption, but nec- 
essarily connected with a spin flip of the electron [^]. 

For T > we see in Fig. [l| that with increasing tempera- 
ture for the spin- J, electron spectral weight is transferred from 
the high-energetic polaron peak to the low-energetic scatter- 
ing peak. On the other hand, for the spin-f electron we have 
the effect that for finite temperatures an additional peak rises 
at the high energetic side of the spectrum. This rise with in- 
creasing temperature is fueled by the loss of spectral weight 
on the low-energetic peak. The high-energetic peak simply 
represents the ability of the spin-t electron to exchange its 
spin with a not perfectly aligned local-moment system, T > 
|[r^]. As a result of the shifts of spectral weight occuring for 
both, the spin-t ^nd spin-| spectra, the spectra for the two spin 
directions approach each other with increasing temperature. 
In the limiting case of T Tq the system has eventually lost 
its ability to distinguish between the two possible spin direc- 
tions because of the loss of magnetization of the underlying 
local-moment system. Hence, for T = Tc the densities of 
states of the spin-J, and the spin-t electron are equal. Another 
feature which can be seen in Fig. |l] is that while the positions 
of the quasiparticle subbands stays pretty much the same we 
observe a narrowing of the bands with increasing temperature. 



This can especially well be seen for the case of temperature 
evolvement of spectrum of the spin-| electron for the case of 
J = 0.2. Whereas the scattering band and the polaron band 
are, for the case of T = 0, still merged, in the case of T = Tc 
both bands are clearly separated. 

In this paper we are interested in surface states and their 
temperature dependent behaviour. Surface states occur in the 
spectral density at energies different from the bulk energies 
and are localized in the vicinity of the surface of a crystal. 
The theory presented above is applied to a s.c. film consist- 
ing of n layers oriented parallel to the (lOO)-surface as drawn 
schematically in Fig. ^ 

The electron hopping in Eq. is restricted to nearest 
neighbours. 



T 



a(3 



(17) 



where A stands for the nearest neighbours within the same 
plane, A - (0,1), (0,1), (1,0), (1,0), and d^f = 5^p5,,. 
T°'^ is the hopping between the adjacent layers a and f3 and 
T"" is the hopping between nearest neighbours within the 
layer a. In order to study surface states we vary the electron 
hopping within the vicinity of the surface. 
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according to Fig. ^, and with 

Til = eiiT, 



\ 
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T T T± 
T^ Tn J 



(18a) 



Ti_ = ti_T. 



(18b) 



Here, e|| and ej^ are considered as model parameters. In re- 
ality the variation of the hopping integrals in the vicinity of 
the surface may be caused e.g. by a relaxation of the in- 
terlayer distance. According to the scaling law T ^ r^^ 
for the t/-electrons 124| a relatively small top-layer relaxation 
Ar/?- may result in a strong change of the hopping integral T. 
Thus e.g. a relaxation of the Gd(OOOl) surface layer of 3-6% 
(cf. and references therein) would yield a modification of 
the hopping integrals of of up to 30%. 

In a previous paper we have been dealing with surface states 
for the exactly solvable case of a single electron in an oth- 
erwise empty conduction band and a ferromagnetic ally satu- 
rated/spin system, T = [|2^]. For this special case we have 
shown that modifying the hopping in the vicinity of the sur- 
face according to Eqs. (18) leads to the appearance of surface 
states in the local spectral density S^" . Modifying the hop- 
ping within the surface layer by more than 25%, i.e. e|| ^ | 
or e|| ^ |, while keeping all the other hopping integrals un- 
changed results in a single surface state at the lower or the 
upper edge of the bulk band [|26|J27|]. This surface state first 
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emerges at the f - and at the M-point from the bulk band and 
from there spreads for larger modifications of e|| to the rest of 
the Brillouin zone. 

On the other hand, when the hopping within the first layer 
remains constant, e|| = 1, but the hopping between the first 
and the second layer is significantly increased, e±_ ^ \/2, then 
two surface states split off one on each side of the bulk band. 
In this case the emergence of the surface states from the bulk 
band is k-independent. Both types of surface states for the 
special case of T = can be observed at the single bulk 
band and on the high-energetic polaron band for the case of 
the spin-t electron and the spin- J, electron, respectively. 

Here, we are interested in the possible variations of the 
T = surface states when going to finite temperatures. 
Figs, ^to H show the temperature dependence of surface states 
for the different possible variations of the hopping in the vicin- 
ity of the surface. All the calculations for Figs. || to ^ have 
been performed for a 20-layer s.c. film cut parallel to the 
(lOO)-plane of the s.c. crystal. For much thinner films than 
20 layers the calculation of surface states becomes meaning- 
less since there is no real bulk-like environment in the center 
of the film to compare the electronic states at the surface to. 
So it is desirable to calculate thicker films for the discussion 
of surface states. The choosen thickness of our model film is 
basically a compromise between computational accuracy and 
computational time. The parameters for the uniform hopping 
according to Eqs. (jl8|) and for the exchange interaction are 
T = —0.1 and J — 0.1, respectively. 

For our calculations we have employed a modification of 
the hopping in the surface layer and between the surface layer 
and the adjacent layer of 50% (e^^ 0.5, 1.5) and 100% 



(ej^ = 2), respectively (cf. Eq. (18b)). Such a drastic mod- 
ification of the hopping integrals in the vicinity of the surface 
is rather unlikely to occur in reality. However, as has been 
shown in [|6|| the actual peak position of a surface states de- 
pends only weakly on the variation of e|| and ej_, respectively. 
The selected strong variation of the hopping parameters Tn 
and T± give rise to pronounced surface states which enable 
us to more clearly see the qualitative behaviour of the surface 
states as a function of temperature. 

In Fig. ^we have the case that the hopping within the sur- 
face layer is enhanced by 50%, leading to the existence of a 
surface state on the outer edge of the bulk dispersion. This 
surface state can most clearly be seen at the F-point and the 
M-point in the two-dimensional Brillouin zone. In Fig. ^ the 
spectral density of the first layer, S'^^(i?), for both of these 
points is displayed as a function of energy. If the tempera- 
ture is increased, we see that the position of the spin-| and of 
the spin-| surface states approach each other in a Stoner-like 
fashion until both peaks are equivalent for T = Tc- The os- 
cillations in the spectral densities, which can be seen in Fig. H 
around -0.6eV, 0.3eV, and 0.7eV are due to the finite thickness 
of our model film, as are the respective oscillations in Figs. ^ 
and I 

Also in Fig. ||, but upside down, the spectral density of one 



of the central layers of the 20-layer film, S'^°^''(i?), can be 
seen again for the f-point and the M-point and both spin di- 
rections indicating that the positions of the surface states vis- 
ible in Sl.]^{E) lie outside of the bulk spectrum of the crystal 
for all temperatures. The same is valid for the spectra dis- 
played in Figs. ^ and || In these figures, however, the local 
spectral densities of the central layers have been omitted for 
clarity. 

In Fig. ^ the temperature dependence of surface states is 
documented for the case where the hopping within the first 
layer is reduced by 50%, e|| = 0.5 (e_L = 1). In this case we 
observe a spin-mixing behaviour where the positions of the 
spin-t and of the spin-| surface states stay the same when the 
temperature is risen but spectral weight is being transferred 
between the different peaks. This results in equal populations 
of the spin-], and the spin-f peaks at T = Tc- 

To round up the picture we see in Fig. || the case where 
the hopping between the first and the second layer is modi- 
fied, e± = 2, while the hopping within the first layer remains 
equal to the uniform hopping within the film, e|| = 1. Here 
we have for T = two surface states, one on each side of 
the bulk spectrum. When the temperature is switched on the 
surface states on the outer side of the bulk dispersion behave 
Stoner-like, while the surface states on the inner side of the 
bulk dispersion exhibit a spin-mixing behaviour 

Apparently, our model is able to reproduce a Stoner like 
collapse of the spin- J, and the spin-t peak positions for Tc as 
well as a spin-mixing behaviour In the spin-mixing case there 
are two peaks which both have a majority- and a minority-spin 
contribution. When the temperature is increased, the spec- 
tral weight of these contributions is altered until for T = Tc 
for each peak the spin-] and the spin-f contributions have the 
same spectral weights. 

In being able to reproduce both Stoner-like and spin-mixing 
behaviour, depending on the variation of the hopping and the 
position in the Brillouin zone, our model calculations are in 
harmony with more recent (inverse) photoemission studies on 
Gd(OOOl) films which abandoned the grasp of earlier works 
that the temperature dependent behaviour of the Gd(000 1 ) sur- 
face state has to be either Stoner-like or of spin-mixing type 
[ [rT| , [l2| ]. Especially, it has been shown here that for certain 
parameters it is possible to observe both kinds of behaviour at 
the same time (see Fig. ||). This feature of our model calcula- 
tion seems to be in strong agreement with a scenario proposed 
by Donath and Gubanka [nin. 
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FIG. 1. Density of states pa{E) of a s.c.-(lOO) double layer for 
different s-f interactions J and different temperatures T (in units of 
Tc)- The dotted lines correspond to the case of vanishing i-/interac- 
tion, J — 0, where there is no distinction between spin-| and spin- J, 
electron, P]{E) — pi{E). The spectra for the increasing values of 
the i-/ interaction, J = 0.1, 0.2, 0.3, 0.4, can be associated by the 
increasing distances of their peaks from the J — 0-spectrum. 
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FIG. 2. Model of an n-layer film with s.c. structure. The nearest 

neighbour hopping is assumed to be Ty within the surface layer, T± 

between the surface layer and the layer nearest to the surface layer, 

and T within and between inner layers. 
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J=0.1, e„=1.5, e,=1 
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FIG. 3. Local density of states, Sii{E) of the first layer of a 
20-layer s.c.-(100)-film at the f -point and the M-point for both spin 
directons and different temperatures, T = 0, 0.7, 0.9, 1 (in units 
of Tc) and modified hopping within the first layer, e]\ = 1.5, while 
ej_ = 1. Upside down but in the same sacale, the local spectral den- 
sity of a central layer, Si^l^^\E), is displayed for the same k-points 
and spin directions. 



FIG. 4. Local density of states, Si]^{E) of the first layer of a 
20-layer s.c.-(100)-film at the F-point and the M-point for both spin 
directons and different temperatures, T = 0, 0.7, 0.9, 1 (in units 
of Tc) and modified hopping within the first layer, e\\ = 0.5, while 
e± = 1. 
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J=0.1, e„=1, e,=2 
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FIG. 5. Same as Fig. ^for modified hopping between the first and 
the second layer, e± = 2, while eii = 1. 
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